Hole Transitions in Multiply-Charged Ions for Precision Laser Spectroscopy and 

Searching for a- variation 
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We consider transitions of electron holes (vacancies in otherwise filled shells of atomic systems) 
in multiply-charged ions that, due to level-crossing of the holes, have frequencies within the range 
of optical atomic clocks. Strong El transitions provide new options for laser-cooling and trap- 
ping, while narrow transitions can be used for high-precision spectroscopy and tests of fundamental 
physics. We show that hole transitions can have extremely high sensitivity to a- variation due to 
large relativistic corrections, and propose candidate transitions that have much larger a-sensitivities 
than any previously seen in atomic systems. 



I. INTRODUCTION 

Recent evidence from quasar absorption spectra sug- 
gests that there may be a spatial gradient in values of 
the fine-structure constant, a = e 2 /hc, across cosmolog- 
ical distance scales [l|. The data come from some ~ 300 
different absorbers using two different telescopes, Keck 
and the Very Large Telescope, which together provide 
coverage of the whole sky. All existing astronomical data 
measuring a-variation are consistent with the existence 
of a constant spatial gradient in a (dipole model) @. 

The existence of a spatial dipole in a could be con- 
firmed using terrestrial clocks [3j . The motion of the Sun 
relative to the measured dipole (that is, towards a region 
of the Universe with larger a) implies a time- variation of 
around a/a ~ 10 -18 — 10 _19 yr _1 . The annual motion of 
the Earth modulates the signal, giving an additional shift 
5a /a — 1.4 x 10~ 20 coswi, where to refers to the frequency 
of the annual orbit. To measure this shift, one must com- 
pare two clocks with different sensitivities to a over the 
course of several years [ij . The best current limit on ter- 
restrial time- variation, a/a = (—1.6 ± 2.3) x 10 _17 yr _1 , 
comes from comparison of Hg + and Al + optical clocks 
over the course of a year Q . Here the Hg + transition is 
strongly dependent on a, while the Al + clock is relatively 
insensitive H, || . 

To confirm the dipole model of a- variation in the labo- 
ratory requires finding systems where the sensitivity to a 
is enhanced. The sensitivity is indicated by q-coefficients, 
usually defined by 



U) = uiq + qx 



(1) 



where x — (a/ag) 2 — 1 « 2(a — ao)/ao is the change 
in a over time from its current value ao, which leads 
to a change in the frequency ui. Here the atomic unit 
of energy, which is cancelled out in any measured ratio 
of frequencies, is assumed to be constant. The depen- 
dence of u) on a is due to relativistic corrections. Poten- 
tial clocks that have large q-values include optical tran- 
sitions in Yb+ and Th 3+ Q, as well as the thorium 
"nuclear clock" [9|, which would make use of the 7.5 eV 
nuclear transition in 229 Th to produce a clock with a q- 
value many orders-of-magnitude larger than the optical 



Hg + clock [T(l|Tlj]. Recently, we showed that trapped, 
highly-charged ions could provide optical atomic clocks 
with much larger (/-values than would be possible using 
near- neutral ions [12j . Clocks using suggested transitions 
in Sm 14+ , for example, would benefit from several sep- 
arate contributions to a-sensitivity: high nuclear charge 
Z, high ionization stage, and significant differences in the 
configuration composition of the states involved. 

In this paper we demonstrate that using holes can dra- 
matically increase sensitivity to a- variation, while retain- 
ing all of the other sources of enhancement. We show 
that it is possible to find hole transitions in highly ion- 
ized atomic systems that are within the range of optical 
lasers, and discuss the sources of the large q values. Two 
promising systems are identified along the ionization se- 
quences of iridium and tungsten. The Ir 16+ and Ir 17+ 
ions are shown to have q-values much larger than any 
that have previously been seen in atomic systems. 



II. THEORY 

In Refs. 0, UH it was shown that for a single electron 
above closed shells, the relativistic shift is 
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where /„ = Z 2 /2v 2 is the ionisation energy of the orbital 
(atomic units h = e = m e = 1). Here Z a is the effective 
charge that the external electron "sees" and v is the ef- 
fective principal quantum number. It was shown in [l2[ 
that, for a valence orbital in an ion with large enough 
Z, q ~ /„ as v — > n, the principal quantum number. 
The ratio g n orm = —q/(Z 2 a 2 I n ) was computed for the 
5s valence orbital along the Ag isoelectronic sequence; 
we reproduce the results (Fig. 1. of [12j) as the filled 
circles in Fig. [1] below. 

Consider how might be modified for a hole. For 
our purposes, the hole state has the same quantum num- 
bers as an electron in the closed shells. The hole state 
will have an enhanced sensitivity to a- variation since the 
ionization potential is larger when there are more elec- 
trons in the outer shell, and is maximal for electrons in 
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FIG. 1. Ratio g n orm = —q/(Z 2 a 2 I n ) for the 5s orbital, cal- 
culated using the Dirac-Fock theory. Filled circles: g n orm for 
the 5s valence orbital in the Ag isoelectronic sequence. Open 
squares: q a0 xm f° r a 5s orbital from the configuration 5s 2 4/ fc 
where the ionization degree is always Zi on = 12. 



a closed shell. Indeed, Z a for an electron in closed shells 
would not be well described by Z a = Z- lon + 1 as it was 
previously; rather such an electron would spend approx- 
imately half its time closer to the nucleus than the other 
electrons in the same shell, and this leads to a larger rela- 
tivistic shift. In such a case q ~ if/ 2 will better describe 
the system. 

The open squares in Fig. [1] show the ratio q norul for 
the 5s orbital along the sequence 5s 2 4/ fc with constant 
Zi on , computed in the Dirac-Fock approximation. The 
sequence starts at 5s 2 (k = 0) at Z = 60, and ends at 
5s 2 4/ 14 at Z — 74. For all ions in the sequence Z\ OTy = 12. 
Along this sequence the 4/ orbital energies lie above the 
5s orbital energies by less than a few percent - in this 
sense the orbitals are in the same "shell" . 

Fig. [T] shows that there is indeed an extra contribu- 
tion to the relativistic shift when the external orbital has 
many electrons in it: q is increasing much faster than I n . 
In the Hartree-Fock approximation for filled subshells, 
the ionization energy (to create a hole) is equal to the 
orbital energy. We see that a single hole in an other- 
wise filled external shell will have the highest sensitivity 
to variation of a possible given a particular Z and Z a . 
It is for this reason that we consider hole transitions as 
candidates for optical clocks. 

To exploit the enhanced sensitivity to ev-variation af- 
forded by hole transitions, we must find examples that 
have high Z, are highly ionized, and lie within the optical 
range. This is not trivial, since the difference in ioniza- 
tion potential between different shells increases as ~ Z 2 , 
which rapidly takes any transition between the shells out- 
side the range of lasers as Z a increases. To combat this 
tendency we must find particular examples near crossing 
points between the shells. The Coulomb crossing hap- 
pens because in neutral atoms shells with larger angu- 
lar momentum, I, can have energies significantly above 
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FIG. 2. Dirac-Fock energies of 4/7/2 (solid) and 5s (dashed) 
orbitals in the Sm isoelectronic sequence, N = 62, with closed 
external shells 5s 2 4/ 14 . 



shells with smaller momentum but larger principal quan- 
tum number, n (for example, in neutral thorium the 5/ 
orbital is above 7s). On the other hand, in the hydro- 
genlike limit all shells with the same n are nearly degen- 
erate, regardless of angular momentum (e.g. E§f = E§ s ). 
In this limit, a higher principal quantum number is nec- 
essarily associated with a larger energy. Therefore as Z 
increases between the two limits there must be a cross- 
ing point where the two shells with different n and I have 
similar energies. Near such a point it may be possible to 
find hole-transitions within the range of lasers. 



III. IRIDIUM 

In this paper we identify two separate crossing points 
occurring in high-Z ions. The first happens between the 
4/ 14 and 5s 2 shells in Sm-like ions (number of electrons 
N = 62) with Z > 70. In Fig. [5] we show the Dirac-Fock 
energies of these closed, outermost shells as a function of 
Z. The crossing is seen to occur around Z = 77 (irid- 
ium). In the Dirac-Fock calculation a single hole has 
energy given by the ionization potential, —E. While the 
ionization potentials of both shells are around 380 eV, the 
transition energy between configurations of Ir 16+ having 
the hole in the 4/ level (5s 2 4/ 13 ) and the hole in the 5s 
level (5s 4/ 14 ) will be very much less than this. Config- 
uration interaction calculations are presented in Table HI 
Note that for Z below 70 the 5p shell can begin to fill 
before the 4/ shell. This is not necessarily a bad thing: 
for some Z and some N < 62 there will be crossing be- 
tween partly-filled 5p and 4/ shells that may be exploited 
to find optical transitions, however we have not consid- 
ered these many-electron cases in this work because no 
additional advantage is gained from more complex cases. 

Table Q] and Table [TT] show configuration interaction 
(CI) calculations for low- lying levels in Ir 16+ and Ir 17+ : 
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TABLE I. Energy levels and sensitivity coefficients (q) for 
lr 16 + (cm- 1 ). 



Configuration 


J 


Energy 


<? 




7/2 










5/2 


25 898 


23 652 


4/ 14 5s 2 S 


1/2 


37460 


367 3f 5 



TABLE II. Energy levels and sensitivity coefficients (q) for 
Ir 17+ (cm -1 ). Uncertainties between terms are at the level 
~ 6000 cm -1 , and in particular it is possible that the ground 
state is actually 4/ 14 ^So. However, since this level is 
metastable, it does not qualitatively affect our analysis of this 
system. 



Configuration 


J 


Energy 


q 






Odd states 




4/ 13 5s 3 F° 


4 








4/ 13 5s : 


ipa 


3 


4 838 


2 065 


4/ 13 5s : 


ipo 


2 


26 272 


24183 


4/ 13 5s : 


Ipa 


3 


31492 


25 052 






Even states 




4/ 14 X S 







5 055 


367161 


4/ 12 5s 2 


3 H 


6 


35 285 


-385 367 


4/ 12 5s 2 


? F 


4 


45 214 


-387 086 


4/ 12 5s 2 


3 H 


5 


59 727 


-362 127 


4/ 12 5s 2 


? F 


2 


68 538 


-378 554 


4/ 12 5s 2 


l G 


4 


68 885 


-360 678 


4/ 12 5s 2 


? F 


3 


71 917 


-362 313 


4/ 12 5s 2 


S H 


4 


92 224 


-339 253 


4/ 12 5s 2 


l D 


2 


98 067 


-363 983 


4/ 12 5s 2 


\J 


6 


110 065 


-364 732 


4/ 12 5s 2 


3p 





110 717 


-372 570 


4/ 12 5s 2 


3p 


f 


116 372 


-362 937 



the one-hole and two-hole cases, respectively. Details of 
the CI calculations have been presented elsewhere (see, 
e.g. (flUH). The q- values are calculated by repeating the 
entire calculation for x — —0.01 and 0.01 and taking the 
gradient of the frequency, i.e.: 

U8x — LU-Sx 

2 ox 

with 5x = 0.01. The resulting q values are the largest 
ever seen in an atomic system. Energy levels are accurate 
to 6000 cm -1 , mostly due to two-electron excitations in 
the CI calculation not being fully saturated. However, 
the fine-structure splitting within terms is much more 
accurate than this, and the g-values are very stable and 
certainly correct to within 10%. 

Let us take two particular examples from Ir 17+ , the 
two-hole case. Trapping and cooling of this ion may be 
facilitated by any of the El transitions from the ground 
state 4/ 13 5s 3 F°, such as 4/ 12 5s 2 3 F 4 at ui ~ 45 000 cm -1 . 
While this transition itself has an extremely large q- 
value, and could therefore in principle be used to test 
a-variation, El transitions tend to be too broad for pre- 
cise measurements. On the other hand the transition to 



4/ 12 5s 2 3 Hq with frequency ~ 35 000 cm -1 could be in- 
duced via hyperfine mixing with its fine-structure partner 
3 i?5, and this transition would have a strongly reduced 
linewidth (it may also go as E3/M2). The most abun- 
dant isotope of iridium is 193 Ir, which is stable and has 
nuclear spin / = 3/2. Therefore the hyperfine-induced 
transitions are always available. 

The 4/ 14 x So line is interesting because the a- 
sensitivity acts in the opposite direction (q is large and 
positive). However, when trying to use this line one ex- 
pects to encounter the problem of too little line strength. 
The best option here is to use the transition from the very 
metastable 4/ 14 % to 4/ 13 5s 3 F 2 ° (u> w 21000cm -1 ). 
This is an M2/E3 transition, and at first seems too 
slow for existing lasers. However, due to the hyper- 
fine mixing with J = 1 states such as the 116 000 cm -1 
4/ 12 5s 2 3 Pi level, this transition may also proceed via 
hyperfine-induced El transition. Odd-parity J = 1 
states also exist; hyperfine mixing with these states would 
be dominated by the configuration 4/ 12 5s5p at energies 
~ 500 000 cm -1 ~ 60 eV, which is still far below the con- 
tinuum of such a highly-ionized system. 

The best current limit on terrestrial time-variation of 
a comes from an experiment comparing Hg + and Al + 
clocks [H, with Aq — 57 000cm -1 [f|. By contrast, an 
experiment that compared the two transitions we have 
examined (4/ 14 % to 4/ 13 5s 3 P 2 ° and 4/ 13 5s 3 F 4 ° to 
4/ 12 5s 2 3 Hg) would have a total sensitivity to a- variation 
of Aq ?» 730 000 cm -1 . The largest q ever found in a sin- 
gle atomic transition occurs in 4/ 14 to 4/ 12 5s 2 3 F%, 
however in this case the upper level has an El-decay 
mechanism and may be too broad. 



IV. TUNGSTEN 

The second crossing point we consider happens when 
the filled 4/ 14 and 5p 6 shells cross in Er-like ions (N = 
68). Because of the relatively large splitting between the 
5p 3 / 2 and 5pi/ 2 orbitals (compared to the much smaller 
fine-structure splitting of the 4/ orbitals), this crossing 
point is more spread out, somewhere between Z = 73 
and 75 (see Fig. |3j) . We focus on Z = 74 because tung- 
sten is a well-studied element due to its use in plasma 
diagnostics (see, e.g. [3), and hence may be of higher 
interest to experimenters. CI calculations for the one- 
hole and two-hole systems, W 7+ and W 8+ respectively, 
are presented in Tables IIIII and IIV1 

Previous studies extrapolating from the measured W 6+ 
spectrum [l6[ have noted that the 4/ and hp orbitals 
compete for the ground state in W 7+ [l5j]. We present 
these calculations in Table IIIII but repeat the warnings 
of [15j that the accuracy may be low. The fine-structure 
splittings are more accurate. 

Both W 7+ and W 8+ have several Ml and E2 transi- 
tions that may be suitably narrow to be of use in studies 
of a- variation. For example, the E2 transition from the 
ground state of W 8 +, 4/ 14 5p 4 3 P 2 , to 4/ 12 5p 6 3 F 4 has 
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FIG. 3. Dirac-Fock energies of 4/ (solid), 5s (dashed), and 5p 
(dot-dashed) orbitals in the Er isoelectronic sequence, TV = 
68, with closed external shells 5s 2 5p 6 4/ 14 . The orbital fine- 
structure splitting is indicated, with closed circles for 4/ 7 / 2 
and 5p3/2 and open circles for the 4/ 5 / 2 and 5pi/2 orbitals. 



TABLE IIL Energy levels and sensitivity coefficients (q) rel- 
ative to the ground state for W 7+ (cm -1 ). 



Configuration 


J 


Energy 


<7 






This work 


[15] 




4/"5p t ' ' l F° 


7/2 













5/2 


18199 


17440 


16 462 


4/ 14 5p 5 2 P° 


3/2 


4351 


800 (700) 


87 544 




1/2 


93 908 


87 900 (700) 


200 269 



V. CONCLUSION 

We have shown that level-crossings between filled shells 
provide an opportunity to find hole-transitions within 
laser range in highly-charged ions. These hole transitions 
can have hugely enhanced sensitivity to a- variation when 
compared with transitions used today. As atomic spec- 
troscopy in electron-beam ion trap s continues to improve 
(see, e.g. [l7l4l"9j and review [2fjj ]) we expect that hole- 
transitions in multiply-charged ions may become good 
systems for searching for time-variation of a. We have 
found candidates in two different systems that may be of 
interest, and suggest that others may exist along other 
isoelectronic sequences. Level crossings also allow us to 
identify optical-range El transitions such as those we 
have found in Ir + that may be of importance in laser- 

TABLE IV. Energy levels and sensitivity coefficients (q) for 
W 8+ (cm -1 ). Ener gies between terms are uncertain at the 
level ~ 6000 cm" 1 . 



4/ i4 5p 
4/ 13 5p ; 
4/ 13 5p s 3 
4/ 1 V 
4/ 13 5p 5 3 
4/ 13 5p 5 
4/ 13 5p 5 
4/ 14 5p 4 
4/ 13 5p 5 
4/ 13 5p 5 
4/ 12 5p 6 3 
4/ 12 5p 



4f 12 5p 6 
4/ 12 5p 6 



•ation 


J 


Energy 


q 


3 P 


2 










4 


6 075 


-81564 


3 G 


3 


6 357 


-81480 


3 G 


5 


11 122 


-82 880 


'•F 


3 


21905 


-66 489 


"'D 


2 


23 276 


-66 985 


3 F 


2 


28112 


-66124 


l S 





29 810 


12 735 


3 G 


4 


34 884 


-65 896 


"'D 


1 


36 497 


-65 946 


3 H 


6 


56 416 


-179 990 


3 F 


4 


65 008 


-181419 


3 H 


5 


73 188 


-164139 


l G 


4 


80 551 


-162 573 


3 F 


2 


82 424 


-177137 


3 F 


3 


83 315 


-164072 



energy uj ~ 65 000 cm -1 and q = —181000 cm" 1 . This 
transition may be compared to the 4/ 14 5p 4 1 S'o transi- 
tion, or to any line with a positive or small q value in 
another ion, to give a sensitivity to a-variation which 
would still be the largest ever utilised. By comparing in- 
stead to the 4/ 13 5p 6 2 F° /2 to 4/ 14 5p 5 2 P° /2 E2 transition 

in W 7+ , one could achieve Aq w 365 000cm" 1 . 



trapping and cooling of these ions. 

ACKNOWLEDGMENTS 

This work was supported by the Australian Research 
Council. Supercomputer time was provided by an award 
under the Merit Allocation Scheme on the NCI National 
Facility at the Australian National University. 



[1] J. K. Webb, J. A. King, M. T. Murphy, V. V. Flambaum, 
R. F. Carswell, and M. B. Bainbridge, "Evidence for 
spatial variation of the fine structure constant," (2010), 
arXiv: 1008.3907. 

[2] J. C. Berengut and V. V . Flambaum, J. Phys.: Conf. 
Ser.. 264. 012010 (2011). larXiv:1009. 3693 



[3] J. C. Berengut and V. V. Flambaum, "Manifestations of 
a spatial variation of fundamental constants on atomic 
clocks, oklo, meteorites, and cosmological phenomena," 
(2010), arXiv:1008.3957. 

[4] V. A. Dzuba, V. V. Flambaum, and J. K. Webb, Phys. 
Rev. Lett., 82, 888 (1999). 



5 



[5] T. Rosenband, D. B. Hume, P. O. Schmidt, C. W. Chou, 
A. Brusch, L. Lorini, W. H. Oskay, R. E. Drullinger, 
T. M. Fortier, J. E. Stalnaker, S. A. Diddams, W. C. 
Swarm, N. R. Newbury, W. M. Itano, D. J. Wineland, 
and J. C. Bergquist, Science, 319, 1808 (2008). 

[6] V. A. Dzuba, V. V. Flambaum, and J. K. Webb, Phys. 
Rev. A, 59, 230 (1999). 

[7] S. G. Porsev, V. V. Flambaum, and J. R. Torgerson, 
Phys. Rev. A, 80, 042503 (2009). 

[8] V. V. Flambaum and S. G. Porsev, Phys. Rev. A, 80, 
064502 (2009). 

[9] E. Peik and Chr. Tamm, Europhys. Lett., 61, 181 (2003). 
[10] V. V. Flambaum, Phys. Rev. Lett., 97, 092502 (2006). 
[11] J. C. Bereng ut. V. A. Dzuba. V. V. Flambaum. and 
S. G. Porsev, |p"hys. Rev. Lett., 102, 210801 (2009) . 



[12] J. C. Berengut, V. A. Dzuba, and V. V. Flambaum, 

Phys. Rev. Lett., 105, 120801 (2010). 
[13] J. C. Berengut, V. A. Dzuba, V. V. Flambaum, and 

M. V. Marchenko, Phys. Rev. A, 70, 064101 (2004). 
[14] J. Clementson, P. Beiersdorfer, E. W. Magee, H. S. 

McLean, and R. D. Wood, J. Phys. B, 43, 144009 (2010). 
[15] A. E. Kramida and T. Shirai, At. Data Nucl. Data Tables, 

95, 305 (2009). 

[16] J. Sugar and V. Kaufman, Phys. Rev. A, 12, 994 (1975). 

[17] I. Draganic, J. R. Crespo Lopez-Urrutia, R. DuBois, 
S. Fritzsche, V. M. Shabaev, R. S. Orts, I. I. Tupitsyn, 
Y. Zou, and J. Ullrich, Phys. Rev. Lett., 91, 183001 
(2003). 

[18] J. R. Crespo Lopez-Urrutia, Can. J. Phys., 86, 111 
(2008). 

[19] M. Hobein, A. Solders, M. Suhonen, Y. Liu, and 

R. Schuch, Phys. Rev. Lett., 106, 013002 (2011). 
[20] P. Beiersdorfer, Phys. Scr., T134, 014010 (2009). 



